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Abstract: 

In an earher paper, we discuss the "popcorn" model for baryon production in quark and gluon 
jets, and present an improved model (which we call Modified Popcorn Scenarium, MOPS). 
In this paper we give a manual to the MC program based on MOPS, and also discuss the 
application of the model to baryon fragmentation, i.e. fragmentation of strings originally 
contaning a diquark. Model predictions for baryon production in DIS are compared with 
data. 
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1 Introduction 



In ref |T[], a model for baryon production in the hadronization process is presented. The 
model is called "Modified Popcorn Scenarium" (MOPS), and is a revision of the "popcorn" 
mechanism |^ for baryon production in the Lund string fragmentation picture In ref 
the discussion is focused on baryon production in e~^e~ annihilation. However, the model 
can also be applied to fragmentation of strings originally containing a diquark. Physical 
considerations about this extension of the model are discussed in section |], where also some 
results are presented. 

A large part of this paper (section ^ and forward) is a manual to a FORTRAN source code 
file called "MOPS.f", which is an implementation of the the new model and works as a 
complement to JETSET> 7.4 The program is intended for those who already have some 
experience of the JETSET Monte Carlo. 



2 Baryon Production 

In the "popcorn" model p[ a diquark is produced in a stepwise manner. In the string it 
is possible to produce a qq pair of "wrong" colour, which does not split the string in two 
singlet systems. If instead the original and new quark(antiquark) form an antitriplet (triplet) 
state, the colour field between the produced quark and antiquark will correspond to a triplet 
colour field with the same strength as the original field. Then no net force is acting on the 
produced quark and antiquark and in accordance with the uncertainty principle they can 
move freely for a time inversely proportional to their energy. If the string breaks within 
the colour fiuctuation region, an effective diquark-antidiquark pair is produced. If several 
break-ups occur, there will be one (or several) intermediate mesons between the baryon and 
antibaryon. 

As no net force is acting on the colour fluctuation quarks they appear as free particles. 
Therefore it was in ref expected that the probability to flnd such a qq pair with transverse 
mass /i_L at a separation d in space is proportional to {Apid, fi±)f . The suppression factor 
for producing effective diquarks and intermediate mesons of a total invariant mass M± in a 
string with energy density k, is thus estimated by 

|A^(rf,/ix)r~exp(-2/ixM^/fi:). (1) 

When implementing the model to a MC, the suppression factor Eq (|l|) was originally rep- 
resented by a set of parameters, governing production rates of different diquarks and the 
probability for one intermediate meson. The small probability for several "popcorn" mesons 
was neglected. In ref the model was reformulated in such a way that Eq ([1|) could be 
used explicitly in the MC. The two parameters 

(3g = 2 (/i^q) /k, or (3^ and A(3 = (3, - /5„ (2) 
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Figure 1: A string with an original diquark can be assumed to have similar properties as a 
string remainder after an effective diquark break-up. 

then govern both the diquark and the intermediate meson production. The result is a reduced 
number of tunable parameters. 

We have assumed that the colour fluctuation quarks appear as free, massive particles. This 
implies that they cannot move along light-cones. Instead there is an area of possible start- 
ing points for the colour fluctuation, which is essentially given by the proper time of the 
production vertex squared. To get the total weight for a colour fluctuation, we must sum 
over all possible starting points. The result is a relative suppression of diquark vertices 
at early times, and consequently a softening of baryon momentum spectra, consistent with 
experimental data. 



3 Diquark Jet Fragmentation 

An algorithm for baron production in a qq string can be applied to fragmentation of a string 
originally containing a diquark. This string is then assumed to behave like the remainder of 
a qq string after an effective diquark vertex, as symbolically shown in Fig |l|. 

This identification must be made with some care. The general statement that a remainder 
string behaves like an original string is based on the fact that two string systems separated 
by a vertex are causally disconnected. In the case of diquark vertices, this is not obviously 
correct. Due to the popcorn mechanism, the string segments A and B in Fig ^ share a 
common set of possible histories. Choosing among those when fragmenting A will influence 
the probability for different final states in B. 

In the specific case of the model presented in the discussion about different production 
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Figure 2: Two string segments A and B, causally connected by a common colour fluctuation 
quark pair. 

points for the colour fluctuation pair, giving rise to the suppression of break-ups at small 
proper times (small-F suppression, where F = /t^r^) is not applicable to fragmentation of 
original diquarks. Furthermore, in the process (original diquark —>■ baryon+q) the decuplet 
suppression should not apply at full strength. This suppression factor stems from the nor- 
malization of the overlapping q and q wavefunctions in a newly produced qq pair ^, but 
in the process considered here, two out of three valence quarks already exist as an initial 
condition of the string. 

Finally the suppression factor for popcorn mesons (which is exponential in the invariant mass 
of the meson system) is derived from the assumption of colour fluctuations violating energy 
conservation and thus being suppressed by the Heisenberg uncertainty principle. When 
splitting an original diquark into two more independent quarks, the same kind of energy 
shift does not obviously emerge. We could still expect large separations of the diquark 
constituents to be suppressed, since the two quarks originate from a common diquark wave 
function. The shape of this suppression is beyond the scope of our model, but for simplicity 
we assume the same kind of exponential suppression as in the "true popcorn" case. However, 
there is little reason to assume that the strength of the suppression must be exactly the 
same in the different situations. Thus the suppression of leading rank meson production 
in a diquark jet is governed by a new /3 parameter, which is independent of the popcorn 
parameters (3^ and Af3. 

To summarize, when implementing the popcorn algorithm to the case of original diquarks, 
we have made the following changes: 

* The F suppression does not apply. 
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* The decuplet suppression does not apply when qq— (rank 1 B)+q, but when qq^ 

M+... + (higher rank B)+q. 

* The P parameter used in the exponential popcorn suppression is independent of /?u and 

A/3, which are fitted to e~^e~ baryon production data. 

To test the model, we have made simulations of 7r+p and K+p collisions at ^/s ~ 20GeV, 
and compared the predictions with the results of ref 0. In the simulations, we have used 
FRITIOF 7.02 and JETSET 7.4 (default and modified). Before presenting the results, we 
want to stress that the low energy-scale involved {y/s ~ 20GeV =^ string energies mostly 
< 10 GeV) makes it dangerous to draw too firm conclusions from iterative Monte Carlo 
simulation programs. 

The results from this preliminary study are briefly (cf Fig |^ and Table |ip: 

* The new model improves fit to A spectra. 

* No version reproduces the "dip" in the spectra. 

* The new model does not reproduce multiplicities as well as the default version. 

* Changing the s/u parameter from 0.3 to 0.2, both versions improve multiplicity predic- 

tions, and now the new model is slightly better than default. 

The last result is of some interest, since HERA data at much higher energies also indicate 
support for s/u~ 0.2 in DIS 0. 
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Table 1: Multiplicities of baryons in K^p collisions at ~ 20 GeV cms energy. Data from ^ 
are compared to simulations with FRITIOF7.02 combined with JETSET default and modified, 
for two different values of the strangeness suppression parameter. N stars indicates that the 
simulated result is more than N standard deviations away from experimental data. 



4 The MOPS Program 

Installation 

All changes of the code are inside the subroutine LUSTRF. The modifications have been 
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Figure 3: Normalized Feynman-x distribution of A and in vr+p collisions at ^20GeV 
cms energy. The meson is moving in the positive direction. Data from ^ are compared 
to simulations with FRITIOF7.02 combined with JETSET default (solid lines) and modified 
with (3 = 0.6 GeV'^ (dashed lines). 



tested with JETSET version 7.403. However, the file MOPS.f does not contain the complete 
JETSET code, but only the modified LUSTRE (and some new subroutines called by it). 
Thus it should be possible to combine the modifications with different sub-versions of JET- 
SET 7.4 and perhaps later versions. To use the modified program, remove the subroutine 
LUSTRE from a copy of the JETSET code. Create object files both of this copy of JETSET 
and of MOPS.f, with the same options you normally use when compiling JETSET, and link 
your main program together with both these object files. The modified program should then 
work exactly as JETSET > 7.4 with all default options available and some new options 
added. 



Two precautions must be made: 



* The new code uses some of the available positions in the LUDATl commonblock. 

// anyone of these has been occupied in your version of JETSET it might be impossible 
to use the modifications. 

* Positions in a commonblock that are unused by JETSET default are given the value 0. 

Thus the new switches and variables must be explicitly initiated by the user. 



Uncovered physical aspects 

This is not an official change of the LUSTRE code. It is not prepared for every possible 
string configuration. Three known shortcomings are: 



Present implementation seriously underestimates baryon production in closed strings (e.g. 
T decay or closed strings from colour reconnection). 
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* Tensor meson production inside a colour fluctuation is not implemented. 

* Modification not available for "junction strings", coming from multiple interactions. 

(This is however a feature not fully implemented in JETSET at present.) 

There are no immediate plans of extending the implementation to cover the points above, 
unless feed-back shows a need to do so. 

Diqueirk probability shift 

In the way the new model is implemented, some of the produced baryons are rejected due 
to r suppression. To get a proper amount of final baryons, the program therefore needs 
a larger -P(qq)/-P(q) value as input. This could of course be achieved simply by changing 
PARJ(l). However, the default flavour selection subroutine, LUKFDI, is sometimes called 
by LUDECY (which handles particle decays). Since the new model is not implemented in 
LUDECY, a drastic enhancement of PARJ(l) could lead to an overestimate of heavy hadrons 
decaying into baryon-antibaryon pairs. 

Instead of including a modified LUDECY in MOPS.f, the shift of -P(qq) / -P(q) is implemented 
as a new parameter, giving a multiplicative factor to PARJ(l) when using the F suppres- 
sion algorithm. This new parameter has however no physical interpretation, and could be 
excluded by implementing modifications to LUDECY. 



5 Subroutines 

Some new subroutines are called inside LUSTRF, but the user should not need to call them 
explicitly. Here is however a short description. 

SUBROUTINE LUSTRF(IP) : fragments a string as default LUSTRF, with options 
added to use F suppression of baryons (MSTJ(12)> 4) and new flavour algorithm 
(MSTJ(12)=5). 

SUBROUTINE POKFDI(KFl,KFQPOP,KF3,KF,PT2,PI5) : generates a new fla- 
vour and combines it with an existing flavour to give a hadron, i.e. the same as 
LUKFDI. POKFDI handles incoming diquarks according to the new popcorn algo- 
rithm, switched on by MSTJ(12)=5. 

KFLl : (input) incoming flavour. 

KFQPOP : (both) flavour of present colour fluctuation quark. When a popcorn 

system is started (q ^ B + qq), KFQPOP is initiated in POKFDI. If KFQPOP 
is a diquark flavour, the subroutine will instead form final baryon, build up by 
KFl and KFQPOP. 

KFL3 : (output) produced new flavour. Is if KFQPOP is a diquark. 

KF : (output) produced hadron. 
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PT2 : (input) p\ of hadron to be produced. 

PI5 : (output) Transverse mass squared of produced hadron. 

SUBROUTINE DIQINI : initiates weights for diquark flavours, baryons and popcorn 
mesons in the MSTJ(12)=5 popcorn algorithm. Is called if MSTJ(5)=0 and sets 
MSTJ(5) to 1. // changing any parameter or switch governing these weights during 
a run, MSTJ(5) must he reset to or DlQlNI called explicitly. 

SUBROUTINE PONMES(IWT) : generates number of mesons to be produced in a 
popcorn system. If the total system is not accepted, a completely new one is started, 
and PONMES is called again. 

IWT : (input) Switch used to identify original diquark fragmentation (=1) or ordinary 
popcorn system (=0). 



6 Parameters and Switches 

The modified program uses MSTJ(5), MSTU(121,122) and PARJ(8-10,20,27). Do not com- 
bine the modification with a JETSET version in which any of these variables are occupied! 

The flavour weights depend differently on the implemented parameters in different situa- 
tions, and the calculations needed to extract proper weights from the parameters can be 
rather lengthy. For efficiency, the program starts a simulation by calculating the large but 
limited number of weights. The result is stored in PARF(201-517). Strictly speaking, these 
positions are occupied: In JETSET, the user has an option to define flavour weights in 
PARF(201-1960). However, since it must be considered nonsense to combine the flavour 
modification of this program with "home-made" flavour weights, this commonblock space 
is used, and a combination of the new flavour algorithm and the MSTJ(15)> option is 
explicitly forbidden. 

Note that the default values (denoted by "(D=...)" below) are completely nonsense! This is 
due to the fact that the subroutine BLOCK DATA which defines default values might change 
significantly between different JETSET versions. In order to make it possible to link the 
modification to different JETSET versions, the BLOCK DATA subroutine is not included 
in the file MOPS.f. Thus the default values are not related to the new model. The prelim- 
inary tune presented in 0] was: MSTJ(5)=0, MSTJ(12)=5, PARJ(8)=0.6, PARJ(9)=1.2, 
PARJ(10)=2.09, PARJ(20)=0.6 (from the preliminary results presented in section ^ of this 
paper) and PARJ(27)=0.19. 

MSTU(121) : (Internal) counts popcorn mesons. Each time a popcorn system is started 
in the new algorithm, the number of mesons to be created is stored in MSTU(121) 
which then is decreased by 1 for each meson produced. When MSTU(121)=0 and 
the popcorn system is accepted in spite of suppressions, a baryon is produced and the 
system completed. 



7 



MSTU(122) : (Internal) switch to remember whether a popcorn system started through 
q ^ B + qq production or a baryon remnant diquark at the string end. 

MSTJ(5) : (D=0) switch to calculate different weights related to the parameters of the 
model. DIQINI is called by LUSTRF if MSTJ(5)=0. Is set to 1 in DIQINI. Reset to 
or call DIQINI explicitly if you change any parameters governing flavour or popcorn 
weights during the run. 

MSTJ(12) : (D=2; suggestion 5) choice of baryon production model. 

= 0-3 : see manual, default version. 

= 4 : popcorn model (as =2) with F suppression. 

= 5 : r suppression and new popcorn algorithm, independent of PARJ(3-7,18). In- 
stead depending on PARJ (8-9,27) (Baryon remnant fragmentation also depends 
on PARJ (20)). 

PARJ(8,9) : (D=0.,0.; suggestion 0.6,1.2) give Pu and A/3 in GeV"^ for the popcorn sup- 
pression factor. 

PARJ(IO) : (D=0.; suggestion 2.1) multiphcative factor to PARJ(l)=P(qq)/P(q), used 
if MSTJ(12)>4, compensating for rejected popcorn systems due to P suppression. 
PARJ(IO) merely corresponds to a shift of the PARJ(l) value. With the present 
implementation of particle decays in JETSET, it is however important not to change 
PARJ(l) far from default. This could affect some decay channels of heavy hadrons in 
an undesired way. 

PARJ(20) : (D=0.; suggestion 0.6) gives /? in GeV~^ for the case of fragmenting a string 
spanned by an original diquark. 

PARJ(27) : (D=0.; suggestion 0.19) gives the suppression of baryons in the decuplet. 
PARJ(27) in the new model is exactly the same as PARJ(18) in the default one. 
The reason for using a new spot in the commonblock is again that a change of the 
flavour parameters used in LUKFDI affects decay channels of heavy hadrons. 
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